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LONG-TERM  GOALS 


The  theory  of  weak  turbulence  is  the  cornerstone  of  modern  wave-predicting  models.  It  is  based  on  the 
Hasselmann  equation  for  the  wave  action  spectrum  &)  : 

8N  dcodN_ 

8t  8k  dr  nl  F  ?  ^nl  =  Sin  +  Sdiss  ^ 

s  s  s 

Here  nl  is  the  nonlinear  wave  interaction  term,  in  is  the  air-wave  coupling  term,  and  diss  is 
responsible  for  the  white-capping  dissipation. 

In  spite  of  tremendous  efforts  invested  in  numerical  studies  of  the  Eq.l,  the  progress  in  development  of 

c 

fast  and  accurate  prognostic  models  is  slow.  Just  as  an  illustration,  the  scatter  of  in  in  different 
models  of  exceeds  its  mean  value  up  to  the  factor  of  4  !.  As  a  result,  the  major  verification  tool  for 
wave-prediction  models  —  hind-casting  of  ocean  surface  wave  events  —  consists  in  “tuning-up”  of 

S  S 

coefficients  in  front  of  in  and  diss  terms  to  achieve  reasonable  correspondence  of  numerical 
simulation  results  and  measurable  field  parameters. 


Our  long-tenn  goal  is  improvement  of  surface  wave  -  prediction  models  through: 
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•  Fonnulation  of  new  in  and  diss  terms,  which  can  be  used  in  wave-predicting  models  without 
necessity  of  empirical  tuning-up  to  achieve  adequate  description  o  the  reality 

•  Search  of  the  universal  self-similar  properties  of  Hassellmann  equation  (1),  describing  the 
behavior  of  wave  turbulence  in  terms  of  limited  number  of  intrinsic  physical  quantities,  such 
as:  wave  energy  spectral  flux,  characteristic  wave  frequencies  and  universal  self-similar 
parameters. 


OBJECTIVES 

g 

Our  first  objective  was  to  concentrate  on  the  study  of  diss  term  through  direct  numerical  comparison 
of  surface  wave  turbulence  in  the  frame  of  full  3D  equations  with  Hasselmann  equation,  for  the 
sitiation  when  the  wind  is  absent,  i.e.  pure  swell  setup. 
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Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 
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The  second  objective  was  to  compare  the  comprehensive  experimental  set  of  data,  obtained  over 
almost  last  50  years  of  field  observations  with  the  result  of  analytical  and  numerical  results  of  study  of 
Hassehnann  equation. 

APPROACH 


For  achievement  of  our  objectives  we  used  the  following  two  numerical  tools: 


•  Statistical  Computational  Wave  Tank  (SCWT)  numerically  simulates  non-stationary 
Hassehnann  equation: 


dN  dco  dN 

dt  dk  dx  F  S„[  =  Sjn  +  Sj 


(3) 


The  algorithm  is  updated  Resio-Tracy  code,  which  is  two  orders  of  magnitude  faster  then  the 
original  version.  The  new  version  is  stable.  We  used  the  grid  in  A-plane  containing  36  angular 

U  —  1  0/72  /  SCC 

versus  7 1  frequency  modes.  For  the  wind  velocity  0  it  takes  eight  hours  of 

computational  time  to  model  two  days  evolution  of  the  wind-driven  sea  in  one  spatial  point. 


•  Dynamical  Computational  Wave  Tanks  ( DCWT)  follows  not  only  the  wave  surface  elevation, 
but  also  the  wave  phases  and  is  based  on  the  truncated  expansion  of  the  Hamiltonian  into  the 
series  by  powers  of  the  surface  elevation: 


rjt  =  ki//-(V(r/V  i//))  -  k[i)ki//]  +  k(r/k[i]kif/])  +  ^  A[ij'  ki//]  +  \ k[r/2  Ai//] 
Yt  ~  ~S71  +  0A77 - y[(V i//)2  -(ky/)2\-[ky/~\k[ijky/\-\ijky/\Ay/ 


where  ,^x,  •''*  is  water  surface  elevation  and  Y(x-> >’)  js  velocity  potential  on  the  surface  of  the 
liquid,  &  and  0  are  gravity  acceleration  and  surface  tension  coefficient.  Action  of  the  operator 

k  =  Jkl 


is  defined  as 


dkdk„ 


-2  +k2 

y  is  the  absolute  value  of  the 


A  =  1  y ,  where 

wave  number.  The  simulation  of  DCWT  was  performed  in  real-space  domain  ^~n  x  "^on  the 
grid  5 12x4096  points. 


We  hired  two  contractors  groups  of  applied  mathematicians  and  physical  oceanographers  in  Moscow: 


1 .  Group  in  the  Shirshov  Institute  of  Oceanology,  Moscow,  Russia.  A  head  of  the  group  is  Dr. 
Sergei  Badulin. 

2.  Group  in  the  Landau  Institute  for  Theoretical  Physics,  Moscow,  Russia.  A  head  of  the  group  is 
Dr.  Alexander  Dyachenko. 
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The  following  are  key  individuals  involved  into  research: 


Principals: 

•  Zakharov  V.E.  -  general  supervision 

•  Pushkarev  A.N.-  Hassehnann  and  2D  dynamical  equations 


Contractors:  Group  1: 


$ 

•  Badulin  S.I.  —  time-limited  nl  code  for  Hassehnann  equation 

$ 

•  Geogdjaev  V.V.  —  time-limited  code  for  Hassehnann  equation 

Group  2: 


•  Dyachenko  A.I.-  ID  exact  dynamical  equations 

•  Korotkevich  A.O.  -  2D  dynamical  equations 

•  Prokofiev  A  -  exact  dynamical  equations 


WORK  COMPLETED 


1 .  Numerical  comparison  of  swell  evolution  in  frame  of  two  models:  SCWT  and  DCWT. 

2.  Comparison  of  the  comprehensive  set  of  20  dependencies,  obtained  over  last  50  years  from 
field  observations,  with  analysis  of  Hassehnann  equation  and  numerical  simulation  of  SCWT 
model. 


RESULTS 

The  first  major  result  is  confirmation  of  the  theory  of  weak  turbulence  through  direct  comparison  of 

S  =0 

SCWT  and  DCWT  models.  Comparison  was  perfonned  for  pure  swell  (no  wind  input,  ) 

situation.  Similarity  was  observed  in  spectral  peak  downshift  and  bi-modality  of  the  spectrum  (see 
Fig.2).  Probability  of  distribution  function  of  wave  amplitudes  exhibits  Gaussian  distribution  (see 

g 

Fig.3).  We  compare  SCWT  and  DCWT  models  for  dissipation  terms  diss :  artificial  viscosity,  WAM 
cycle  3,  WAM  cycle  4  and  new  dissipation  terms.  The  result  of  comparison  of  total  energy  behavior  is 
presented  on  Fig.  4.  It  has  been  learned  that  none  of  the  dissipation  terms  except  new  one  give  good 
correspondence  in  total  energy  behavior  in 

Sdiss  =  GrfsCD  ~~ 

rC 

cds  =  i 


time.  The  formula  for  new  dissipation  tenn  is  the  following: 

((|-*)  +  *l)fe)  nk  (5) 

.00  x  10“6,  <5  =  0,  p  =  12 
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[  Fig.  2  Original  spectrum  (first  picture)  down-shifts  in  DCWT  (second picture)  and 

SCWT  ( third  picture)  models] 


1°5  F  ' 

Time  =  3.7  min,  or  67.1  Initial  Wave  Periods 


[  Fig.  3  Probability  distribution  function  of  wave  amplitudes  is  approximated  by  straight 
line,  which  confirms  Gaussian  distribution  law  ] 
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[  Fig.  4  Total  energy  behavior  for  DCWT  and  SCWT  models  with  different  dissipation 
terms.  The  best  correspondence  reached  for  new  dissipation  term] 


The  second  major  result  is  in  obtaining  of  an  asymptotic  weakly  turbulent  relation  for  the  total  energy 
s  and  a  characteristic  wave  frequency  a>* : 


SCO* 


g 


r  col  dz 


=  oc, 


\  1/ 

\/3 


dt 


r 


(6) 
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The  self-similarity  parameter  was  found  in  the  numerical  duration-limited  experiments  and  was 
shown  to  be  naturally  varying  in  a  relatively  narrow  range,  being  dependent  on  the  energy  growth  scale 
only.  The  analytical  and  numerical  conclusions  have  been  verified  versus  known  field  dependencies 
for  wave  energy  growth  and  peak  frequency  down-shift.  A  comprehensive  set  of  more  than  20 
dependencies,  obtained  over  almost  50  years  of  field  observations,  was  analyzed.  They  demonstrate 
that  the  weakly  turbulent  law  has  a  general  value  and  describes  the  wave  evolution  well,  apart  from  the 
earliest  and  full  wave  development  stages  when  nonlinear  transfer  competes  with  wave  input  and 
dissipation. 


IMPACT/APPLICATIONS 


Potential  future  impact  for  obtained  results  consists  in: 

•  Using  correct  dissipation  term  (5)  in  operational  wave-prediction  models 

•  Using  relationship  (6)  for  surface  wave  turbulence  parameters  forecasts. 


5 


RELATED  PROJECTS 


US  Army  Corps  of  Engineers  MORPHOS  project.  Our  contribution  consists  in  continuation  of  the 
development  of  new  forms  of  wave-breaking  within  the  context  of  nonlinear  energy  fluxes  in  a  wind  - 
wave  spectrum. 
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